Abstract. The X-ray emission from the young supernova rem- 
nant (SNR) SNR N 103B is measured and analysed using the 
high-resolution cameras and spectrometers on board XMM- 
Newton and Chandra. The spectrum from the entire remnant 
is reproduced very well with three plasma components of kT c 
= 0.55, 0.65, and 3.5 keV, corresponding roughly to line emis- 
sion by the O-K, Fe-L, and Fe-K species, respectively. Nar- 
row band images reveal different morphologies for each com- 
ponent. The kT c = 0.65 keV component, which dominates the 
emission measure (4.5 x 10 65 m -3 ), is in ionisation equilibrium. 
C^) This provides a lower limit of 3000 yrs to the age of the rem- 
<^> nant, which is considerably older than the previously assumed 
OA , age of the remnant (1500 yrs). Based on the measured energy 
?H of the Fe-K feature at 6.5 keV, the hot (3.5 keV) component is 
■ found to be recently shocked (^200 yrs) and still ionising. The 
' high elemental abundances of O and Ne and the low abundance 
of Fe could imply that SNR N 103B originated from a type II 
i—\ " supernova (SN) rather than a type la SN as previously thought. 
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1. Introduction 



2. Observation and Data Reduction 



The young supernova remnant (SNR) N 103B (also known as 
SNR 0509-69) is the fourth brightest X-ray remnant in the 
Large Magellanic Cloud (LMC). It was first identified as a su- 
pernova remnant by Mathewson and Clarke (1983 1 based on 
the non-thermal nature of its radio spectrum and the relative 
emission-line strengths of S II and Ha. Due to its small size 
(3 pc), Hughes et al. ( |1995[ ) estimated the age of SNR N 103B 
to be about 1500 yrs. Recent Chandra data (Lewis et al. 2001) 



show that the remnant is much brighter on the western limb 
than towards the East, this is most likely because of density 
contrasts between East and West. The radio and X-ray bright- 
ness variations are correlated on large s cales, but do not overlap 
on finer spatial scales (Hughes 2001 ). In the analysis of the 
moderate resol ution CCD ASCA spectrum of SNR N 103B, 
Hughes et al. ( |1995| ) found strong emission lines from Si, S, 
Ar, Ca, and Fe, while no emission from O, Ne, or Mg was re - 
quired in order to fit the ASCA spectrum. Hughes et al. ( 1995 ) 
thus concluded that SNR N 103B is the result of a Type la su- 
pernova. However, its proximity to the star cluster NGC 1850 
and the H II region DEM 84 suggest th at SNR N 103B has an 
early-type progenitor (Chu & Kennicutt 1988), providing con- 
tradictory evidence for the nature of the progenitor. 

In this work, we present the X-ray spectra of SNR N 103B 
measured by the XMM-Newton suite of scientific instruments, 
i.e. the Reflection Grating Spectrometers (RGS) and the Euro- 
pean Photon Imaging Cameras (EPIC). Particularly, we exploit 
the high dispersion of the RGS and its unique capability to re- 
solve line emission from extended sources to investigate the na- 
ture of the hot plasma in SNR N 103B in unprecedented detail 
and to shed new light on its progenitor. The large effective area 
of the EPIC cameras at high energies also allows us to study the 
important Fe-K emission. All these are complemented by the 
sharp images produced with the high angular resolution of the 
Advanced CCD Imaging Spectrometer (ACIS) on board Chan- 
dra, which reveal the temperature dependent morphology of 
SNRN103B. 



SNR N 103B was observed by XMM-Newton on 7 August 
2000. The effective exposure times and filters used on the dif- 
ferent instruments are summarised in table [|. The RGS data 
were processed using custom software that is identical in func- 
tion to the RGS branch of the Science Analysis Sy stem (SAS), 
which was described briefly by Rasmussen et al. ( 2001 ). It in- 
cludes pixel by pixel CCD event offset subtraction using me- 
dian diagnostic frames, position dependent CTI correction, and 
correction of dispersion angles imparted by aspect drift. 

Table 1. XMM-Newton observation data of SNR N 103B 



Instrument 


Filter 


Exposure Time 






(ks) 


MOS 1 


Thick 


8.0 


MOS 1 


Medium 


8.0 


MOS 2 


Medium 


23.5 


PN 


Medium 


9.5 


RGS 1 




32.3 


RGS 2 




32.4 
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We chose this analysis path since our response matrix gen- 
erator features support for extended targets, and details in the 
angular distribution of SNR N 103B were included in the re- 
sponse matrices. A single, broad-band archival ROSAT HRI 
image was used, along with mean spacecraft aspect and our 
spectral extraction angular width, to predict the monochro- 
matic response in the readout channel space. This modelling 
approach breaks down as soon as the target's morphology de- 
pends strongly on wavelength. However, it still permits accu- 
rate measurement of total line intensities and significant veloc- 
ity distributions that distort line profiles in the RGS. In the case 
of SNR N 103B, wavelength dependent variations in morphol- 
ogy appear on scales comparable to the XMM-Newton mirror 
point spread function (PSF), and our modelling approach re- 
mains reasonable for the purposes described. However, we do 
still expect to observe small (~20 mA) wavelength shifts for 
some lines due to variations in the remnant's morphology. 
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The raw EPIC data were initially processed with the XMM- 
Newton SAS. This involved the subtraction of hot, dead, or 
flickering pixels, as well as removal of events due to electronic 
noise. The spectrum was extracted from a circular region with 
a radius of 1.5'. Background subtraction for EPIC was done 
using an exposure of the Lockman hole, with similar data se- 
lections and the same extraction region as for the SNR, scaling 
according to the exposure time. Periods of high particle back- 
ground were rejected based upon the 10 — 12 keV count rate 
for the entire field of view. 

In order to support our spectral analysis of the entire rem- 
nant with spatially resolved information on the morphology, we 
also used archival Chandra data of SNR N 103B, which was 
observed with the back-illuminate d ACIS -S3 chip on 4 De- 
cember 1999 (see also Lewis et al. 



2001) 



The cleaned event 
list was taken from the public Chandra archive. Narrow band 
images were extracted with energy selections made according 
to the nominal energy column, which provides a good energy 
estimate of each recorded photon. 

3. Data Analysis 

3.1. XMM-Newton Spectra 

The RGS spectrum of SNR N 103B is shown in figure |l} The 
spectrum is dominated by emission lines from highly ionised 
species of O, Ne, Mg, Si, & Fe. No emission longward of 23 A 
is detected. This can probably be attributed to line-of-sight ab- 
sorption towards the LMC rather than to the absence of N and 
C species, which could have been readily detected by the RGS. 
An important point to note is that we clearly detect O, Ne, and 
Mg lines. The Ovil and Ovill lines are particularly strong. 
These features were not detected in the CCD spectrum of SNR 



N 103B obtained with ASCA (Hughes et al 1995), nor in a 



more recent CCD measurement with Chandra ACIS (Lewis et 



al. [2001! ). In order to observe shorter wavelengths than pro- 
vided by the RGS we also employ the EPIC-MOS CCD spec- 
trum, which is presented in figure |[ Emission lines from He- 
like S, Ar, Ca, as well as an Fe-K blend at 6.5 keV are clearly 
observed. 

The spectral analysis was p erform ed using the SRON 
SPEX package (Kaastra et al. 1996|), which contains the 



MEKAL atomic database (Mewe et al. 1995 ) for thermal emis- 
sion. Both Coronal Ionisation Equilibrium (CIE) and non- 
equilibrium ionisation (NEI) models are available in SPEX. 
Our model for the X-ray emission of SNR N 103B consists of 
several plasma components. For each component, we fit for the 
volume emission measure (n e ntiV), the electron temperature 
(T c ), the elemental abundances, the redshift of the source, and 
the Hydrogenic column density Nn of neutral absorbing gas 
along the line of sight. In NEI models, the ionisation age (n e t) 



(Kaastra & Jansen 1993) is also a free parameter. Here, n e is 



the electron density, ran is the hydrogen density, V is the vol- 
ume of emitting gas and t is the time since the material has 
been shocked and heated to its current temper ature. A distance 
of 51 kpc to the LMC is assumed (Feast et al. 2002 ). 
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Fig.l. The fluxed RGS of SNR N 103B spectrum. The most 
prominent lines are labeled. 
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Fig. 2. The EPIC-MOS of SNR N 103B spectrum. The solid 
line is the best-fit 3 component NEI model. The most prominent 
line blends are labeled. 



In the following, we describe the procedure used to fit the 
total spectrum emitted by the whole remnant. The entire proce- 
dure was carried out for both CIE and NEI scenarios for com- 
parison. We fitted the RGS spectrum first, since it has the most 
details and therefore is most constraining. Two plasma compo- 
nents are sufficient for obtaining a good fit to the RGS spec- 
trum. Assuming that most of the X-ray emission is due to gas 
with more or less similar composition, we coupled the elemen- 
tal abundances in the two components. This also minimises the 
number of free variables in our fit, which makes it more robust. 

It is clear from the EPIC spectrum that SNR N 103B also 
has a hot component that is hard for the RGS to distinguish. 
In order to model this component, we fit the CCD spectrum of 
EPIC with three plasma components. We coupled the elemental 
abundances of O, Ne, Mg, Si, S, Ar & Ca of the third compo- 
nent to those of the two cooler components. The Fe abundance, 
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however, is left free since we cannot fit the Fe-L and Fe-K with 
the same abundance value. The electron temperature, ionisation 
age and abundances of the two cooler components are fixed 
at the values obtained by the RGS fit. We note, however, that 
the hot component is not very sensitive to the other two com- 
ponents even when all three components are fitted simultane- 
ously. In principle, the EPIC spectra could be fitted well with 
only two plasma components, i.e., without the relatively cold 
one. This explains why previous CCD observations did not de- 
tect the low-temperature component, which is dominated by O 
emission. It is only because of the spectral resolution of the 
RGS and its sensitivity to soft X-rays (> 20 A) that we are able 
to resolve the O VII emission and hence need a cool or under- 
ionised component in the model. Finally, for consistency and in 
order to obtain the ultimate fit, we return to the RGS spectrum 
and refit it with three components. The hot component, which 
is constrained by the EPIC, but hardly produces line emission 
in the RGS spectrum is fixed at its EPIC values. The final pa- 
rameters obtained for the three components are presented in 
Table 2 for both the CIE and NEI cases. We obtain a best-fit 
column density of 7V H = (2.5±0.7) x 10 25 mT 2 . We also ob- 
tain a redshift of 340±70 km s -1 , which is consistent with the 

radial velocity of the LMC (i>LMC = ' ml s_1 )- The S 'S" 
nificance of NEI as well as the elemental abundances resulting 
from the fits are discussed below in §4. 



model for the 2p-3s lines of Fe XVII at ~17 A. This is a well 
known problem and has been observed in many coronal sources 
and in the laboratory. Therefore, it probably has nothing to do 
with NEI effects, although in principle NEI could enhance the 
17 A lines with respect to the 15 A lines by means of inner- 



shell ionisation of Fe XVI. Doron & Behar (2002) suggest a 



model that could rectify the standard collisional-excitation line 



powers in MEKAL and in APEC (Smith et al. 2001). They 



show that including resonant excitation as well as recombina- 
tion and ionisation processes (in CIE) from neighbouring ions 
in the models considerably enhances the 2p-3s lines and conse- 
quently produces the observed line ratios. Most of the other dis- 
crepancies in figure |] have local asymmetrical profiles, which 
are a result of a combination of the following: (i) difference be- 
tween the broadband spatial profile used in the response and the 
actual spatial profiles of individual emission lines, (ii) slight er- 
rors in model or data wavelengths, (iii) slight errors in the line 
spread function which does not account for the real asymmet- 
rical profile of the RGS line shapes. The first of these sources 
of errors (i.e. the spatial profile used) should, however, domi- 
nate. There is also a slight overestimation of the continuum (in 
the range 18-23 A) in the models, which is most likely due to 
uncertainties in the background subtraction. These inconsisten- 
cies between data and models are small and do not affect the 
conclusions drawn in this paper. 



Table 2. Spectral fitting results. Components 1 & 2 reflect the 
RGS results, since the RGS is more robust at constraining the 
parameters in this spectral range. Component 3 is required to 
account for the Fe-K emission in the EPIC Data. 



Parameter 


Covered by 


CIE 


NEI 


Component 1: 

n c nHV(10 64 m* 3 ) 
kT c (keV) 
n c t (10 16 m" 3 s) 
redshift (kms -1 ) 


RGS 


6.2±0.8 
0.21±0.05 

340±70 


2.6±0.1 


U.JJ_ 32 
9 -3 + 10 

340±70 


Component 2: 

n c n H V (10 64 m -3 ) 
kT c (keV) 
n c t(10 16 m -3 s) 
redshift (kms -1 ) 


RGS 

+ 
EPIC 


65.1±1.1 
0.65±0.05 

340±70 


65.2±1.2 
0.65±0.05 
>250 
340±70 


Component 3: 

n c n H V (10 64 m -3 ) 
kT c (keV) 
n c t (10 16 m -3 s) 
redshift (kms -1 ) 


EPIC 


4.1 ±0.4 
3.5±0.5 

4100±300 


4.1±0.5 
3.5±0.5 
5.3±0.8 
340±70 


x'/do.f 




1810/902 


1790/899 



The best-fit spectral NEI model for SNR N 103B is com- 
pared with the EPIC-MOS and RGS spectra in figures || and 
^, respectively. The CIE model fits the data nearly as well. 
The overall fit to the data is very good. For the most part, the 
fit residuals are consistent within instrumental calibration un- 
certainties. However, some small discrepancies remain in the 
fit. The most prominent discrepancy is the missing flux in the 



3.2. Chandra Images 

In figure |] we present three narrow -band images obtained with 
the ACIS camera. The images are produced in the energy 
ranges of 0.5-0.7, 0.7-1.0, and 3.0-6.8 keV. The first image 
(0.5-0.7 keV) is essentially an image of O-K emission (O VII 
and O VIII). The second image (0.7-1.0 keV) represents pre- 
dominantly Fe-L emission and the third image (3-6.8 keV) re- 
flects the morphology of Ar-K, Ca-K, and Fe-K lines, as well as 
some continuum. We have also separated this high-energy band 
into its various sub-component images (Ar, Ca, and Fe). How- 
ever, within the statistics available for these images, all their 
morphologies appear to be fairly similar. The three images in 
figure [| correspond roughly to the three plasma components 
described in Table |lj in respective order. For short, we will 
henceforth refer to these three energy bands (and their respec- 
tive plasma components) simply as the O-K, Fe-L, and Fe-K 
bands (components). 

All three images in figure [| show the same asymmetry, 
where the western hemisphere is much brighter than the eastern 
hemisphere. Aside from that, they are quite different in mor- 
phology. The O-K emission is concentrated in the northwest, 
while the Fe-L emission is more evenly distributed along the 
western rim of the remnant. The Fe-K image reveals a bright 
emission spot located at the southwestern rim. We also created 
a continuum subtracted, deconvolved, XMM-Newton Fe-K (i.e 



6.3-6.8 keV) image (see van der Heyden et al. 2002 for details). 
The XMM-Newton Fe-K image (figure ^) also shows a bright 
emission spot located at the southwestern rim, but extends fur- 
ther inward than seen in the Chandra 3-6.8 keV image. How- 
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Fig. 3. The RGS spectrum in the range 1 1-23 A. The solid line represents a best-fit 2 component NEI model (kT c =0.65 & 0.55 
keV) plus an underlying kT c =3.5 keV component as obtained from the EPIC-MOS fits. The overall fit to the data is good. The 
missing flux in the model for the 2p-3s lines of Fe XVII at ~17 A is a known problem which is due missing transitions in the 
plasma code (see text). 




Fig. 4. Chandra images for the approximate energy range 0.5-0.7 ke V (left), 0.7-1 .0 keV (centre) and 3-6.8 keV (right). These 
three energy bands are dominated by emission from O VII & ionOviii (0.5-0.7 keV), Fe-L (0.7-1.0 keV)and Ar, Ca, Fe-K and 
continuum (3-6.8 keV). 
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Fig. 5. The deconvolved XMM-Newton Fe-K image. The Fe-K 
image reveals a bright emission spot towards the southwestern 
rim of the remnant. 




Radius (arcscc) 



Fig. 6. Normalised radial emission profiles for three separate 
photon-energy bands: 0.5-0.7 (dashed curve), 0.7-1.0 (solid 
curve), and 3-6.8 keV (dotted curve). 



ever, while the statistics are of the XMM-Newton image is supe- 
rior the ability to extract more information on the morphology 
of the remnant is hampered by the lower spatial resolution of 
the XMM cameras. 

A more subtle difference between the emission in the three 
different bands is their radial distribution. In order to obtain a 
further insight into the radial distribution, we plotted the nor- 
malised radial profiles in counts per second in the three bands. 
These profiles in counts per pixel are presented in figure ^. 
Each profile is normalised separately to the average count rate 
measured in its energy range within a radius of 17.5". Despite 
the differences in azimuthal distribution, it can be seen that the 
(normalised) brightness of all three components is surprisingly 
similar; peaking at approximately 4", 7", and 10". The Fe-L 
component (dashed curve) lies slightly (0.5"), but consistently 
inside the other two components. Normalising the radial pro- 
files to their peak values does not alter this conclusion. 

4. Discussion 

4.1. Ionisation equilibrium 

As demonstrated by the fitting results in Table ^, the RGS 
spectrum can be described similarly well by either a set of 
NEI or CIE plasma components, giving x 2 /d.o.f. of 1790/902 
and 1810/899, respectively, for the two methods. Moreover, 
both models yield very similar spectra. In terms of emission 
measure, the X-ray spectrum of SNR N 103B is dominated 
by the Fe-L component (2). This emission component is rep- 
resented by an electron temperature of kT c = 0.65 keV in 
both models. The NEI model yields a very high ionisation age 



(n c t > 2.5 x 10 m 6 s) for this component, which is suffi- 
ciently high for the plasma to be in or very close to equilibrium 
(see Me we |2002|). 



For the weaker and cooler O-K component (1), which pro- 
duces the O VII and most of the O VIII emission, the two 
methods actually give very different temperature and emission 
measure values. In CIE conditions, these O ions form at a rel- 
atively low temperature (0.2 keV in the best-fit CIE model). In 
the case of NEI, on the other hand, this component is much 
hotter (0.55 keV), but still ionising (n c t~2.3x 10 16 m~ 3 s). 
Both models reproduce the observed line intensities equally 
well. Note that the very large errors associated with the tem- 
perature (-0.32 keV) and ionisation time (+lel7 m 3 s) in the 
NEI model are due to the fact that the fitting routing finds an- 
other minimum at ~0.2 keV in the case of equilibrium ion- 
isation, as expected. A comparison of the best-fit models for 
the O VII lines in both approaches is shown in figure [7| The 
differences between the NEI and CIE fits can be seen to be 
small, which implies that the data do not decidedly favour one 
scenario over the other. This CIE/NEI ambiguity in the O VII 
lines was recognised before in the analysis of the RGS spec- 



trum of the LMC remnant N132D (Behar et al. 2001 ). It should 



be noted that the ambiguity of NEI/CIE in the O-K component 
(1) does not affect the fits for the other components at all, i.e., 
we get the same parameters for components 2 and 3, regardless 
of whether component 1 is in NEI or CIE. 

The most outstanding evidence for NEI conditions emerges 
from the hot Fe-K component (3), which is responsible for the 
Fe-K feature measured at 6.5 keV, as displayed in figure || 
An attempt to model this feature as a CIE He-like Fe line 
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Fig. 7. Comparison of NEI and CIE fits to the O VII Helium- 
like triplet. The crosses are the RGS data points in the 21.3-22.3 
A range. The solid line is the kT c = 0.55 keV NEI model while 
the dashed line is the kT c — 0.2 keV CIE model. Both models 
provide fits of similar quality, making it difficult to distinguish 
between the models. 



(restframe energy 6.7 keV) requires that it be redshifted by 
~ 4100 km s _1 with respect to the LMC systemic velocity. 
Such a large redshift is not observed for the Ar-K nor the Ca- 
K lines, which arise from the same component as Fe-K and 
overlap spatially. This makes the large redshift interpretation 
unlikely. An alternative and more plausible explanation for the 
spectral position of the Fe-K complex would be a relatively 
low ionisation age of n c t~5.3xl0 16 m~ 3 s for the hot com- 
ponent, as obtained by our NEI fit (see Table [|) without any 
additional redshift. This NEI component predicts a distribution 
of Fe charge states that peaks between Fe XIX and Fe XXIV. 
This component also contributes most of the flux in the Ar and 
Ca lines and provides good fits to these lines. 

4.2. Density, age, and dynamics 

The electron density can be calculated from the emission mea- 
sure (n e riYiV) if we know the volume of the emitting region. 
For simplicity, we generally assume the emitting volume to be 
a thick spherical shell. Based on the emission maps provided in 
figure [| and on the radial profiles shown in figure |[ the shell is 
taken to have a radius of R = 15" and to be AR = R/10 = 1.5" 



thick (Lozinskaya 1992). At a distance of 51 kpc, it gives a 
volume of V = 1.7x 10 51 x/ m 3 , where / is a volume filling 
factor. This is a very crude estimate of the volume, but it is 
sufficient for our order-of-magnitude estimates that follow. 

For the dominant Fe-L component (2), which seems to 
cover most of the western hemisphere (see figure Q), we assume 
a half filled shell or / = 0.5, which gives an electron density of 
n c ~ 25xl0 6 m -3 . Using this density and the ionisation age 
(n c t) from Table 2, we get for the time since the Fe-L gas was 
shocked: t > 3000 yrs. This time is considerably longer than 
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Fig. 8. The EPIC-PN spectrum in the range 5.5-7.5 keV, show- 
ing the Fe-K plot. The solid line is a kT c = 3.5 keV model fit. 
The location of the line centroid (~6.5 keV) indicates that the 
gas is not in ionisation equilibrium. 

the previously assumed age of the remnant (1500 yrs), how- 
ever, the age derived from the ionisation time assumes constant 
density and does not take into account the expansion of the 
shell. 

Since the Fe-K component is concentrated in the south- 
western region of the remnant (see figure Q), we estimate the 
filling factor to be / = 0.25. This gives an electron density of 
~ 8xl0 6 m~ 3 , which in turn yields an age of ~200 yrs. This 
is a relatively short time and is particularly in contrast with 
the derived age of the Fe-L component. For component 1, the 
bulk of the O-K emission is located towards the NW rim of 
the remnant. We again assume a filling factor of / = 0.25. 
This yields electron densities of 7.1 xlO 6 and 10.6xl0 6 m -3 
for the NEI and CIE models, respectively. If the NEI model is 
preferred, its ionisation time would then also imply a time of 
about 100 yrs since shock heating. Note the comparable ages 
of the O-K and Fe-K components. On the other hand, the dif- 
ferent locations (figure [I]) and vastly different temperatures of 
these two components strongly suggest that they may not be 
dynamically connected. 

Indeed, each of the three spectral components can be iden- 
tified spatially. As discussed earlier and as seen in figure |], the 
various spectral components show substantial variation in mor- 
phology. Given our spectral analysis and due to the fact that 
roughly comparable densities are obtained for all of these com- 
ponents, it is most likely that these variations are due to gradi- 
ents in temperature and ionisation ages. The different locations 
and temperatures of the three components suggest that they are 
probably unrelated dynamically. Despite these morphological 
differences the azimuthal average emission profiles show a re- 
markable similarity. Even more surprising is that, compared to 
the Fe L component, the coolest (O Vll/O VIII) and the hottest 
component peak at a larger radius, extend further out, and seem 
to be the most recently shocked components. This would sug- 
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gest that these components are associated with material re- 
cently shocked by the blast wave, whereas the Fe L compo- 
nent may be plasma heated by the reverse shock. However, the 
overabundance of Ar, Ca and Fe K suggests that, at least the hot 
component consists of ejecta material. Note that this confusing 



situation is corroborated by Lewis et al. (2002), who, based on 
imaging spectroscopy with Chandra, show that these elements 
have radially outward increasing abundance gradients. Both the 
abundance gradient and the radial emission profiles of the three 
components do not fit the standard description of young super- 
nova remnants, in which there is an outer, hot, shocked ISM 
component, and a cooler shock ejecta component. 



4.3. Elemental abundances and type of progenitor 

As mentioned previously, SNR N 103B was designated to be 
the result of a Type la SN on the basis of the lack of O, Ne , 
and Mg emission in the ASCA spectrum (Hughes et al. 1995 ). 
The RGS however, is able to resolve prominent O, Ne, and 
Mg emission. Incidentally, the derived abundances are insensi- 
tive to the type of model we use (NEI or CIE). This strongly 
suggests that SNR N 103B could in fact be the result of a 
core collapse SN. To investigate this further, we compare the 
abundance yields as obtained from fits to the RGS spectra to 
those predicted by theory for Type la and Type II nucleosyn- 
thesis models. We also compare these abundances to the av- 
erage LMC values to verify that the observed X-rays do not 
originate from shocked LMC ISM. All LMC abundance values 
are taken from Hughes et al. (1998) except for those of Ar and 
Ca, which are from Russell & Dopita (1992) as these values are 
not supplied by Hughes et al. (1998 ). We use the relative abun- 
dances because they are more consistent among the various 
models than the absolute abundances (see Hwang et al. 2000). 
The calculated abundances for the Type la models, the clas- 
sical W7 and delaye d deto nation (WDD2) model, are taken 
from Nomoto et al. ( 1997 ), while the Type II models for 12 
and 20 M Q zero-age main -sequence progenitors are taken from 
Woosley & Weaver ( 1995 ). The ratio comparisons are given in 
Table |[ The Fe-L abundance is derived from the kT c = 0.65 
keV component, while the Fe-K is from the kT c = 3.5 keV 
component. All abundances are norma lised to Si and relative 
to solar ratios (Anders & Grevese 1989). 



The derived abundance ratios measured for SNR N 103B 
do not compare well with those of the LMC. While some el- 
ements are underabundant by a factor of a few (Ne, Mg, Fe), 
others are grossly overabundant (S, Ar, Ca). The absolute abun- 
dances (not shown) are also larger than those of the LMC, indi- 
cating that the observed emission is by and large due to shocked 
ejecta. Furthermore, the measured abundance ratios compare 
much better with the type II models than with the type la 
models. Particularly, our abundance ratios for O/Si, Ne/Si, and 
Mg/Si are almost an order of magnitude larger than those pre- 
dicted for a type la explosion, while the Fe/Si ratio is about a 
factor of 3 lower than predicted by the W7 model and slightly 
lower than predicted by the WDD2 model. 



Table 3. Element abundance ratios as derived from NEI fits to 
the data compared to the LMC abundances and theoretical yield 
for Type la and Type II (for 12 & 20 M Q progenitors) models. 
The ratios are all relative to solar (Anders & Grevesse 1989| ) 
and normalised to Si. The abundances for the Type la W7 and 



WDD2 models are from Nomoto et al. (1997), while the Type 



II yields are taken from Woosley & Weaver (1995). The LMC 



values are from Hughes et al (1998), except for that Ar and Ca 



which are from Russell & Dopita (1992) 



Abun. 
ratio 


Observed 


Type 1A 


Type II 


LMC 


W7 


WDD2 


12M Q 


2OM 


O/Si 


0.45±0.08 


0.070 


0.019 


0.17 


1.14 


0.61 


Ne/Si 


0.37±0.07 


0.005 


0.001 


0.13 


1.08 


0.93 


Mg/Si 


0.25±0.06 


0.061 


0.019 


0.13 


2.07 


1.03 


S/Si 


1.91 ±0.21 


1.074 


1.167 


1.63 


0.5 


1.16 


Ar/Si 


2.56±0.64 


0.751 


0.935 


2.19 


0.35 


0.54 


Ca/Si 


2.91 ±0.87 


0.891 


1.382 


1.75 


0.42 


0.34 


FeL/Si 


0.52±0.14 


1.88 


0.866 


0.23 


0.30 


1.16 


FeK/Si 


1.24±0.44 


1.88 


0.866 


0.23 


0.30 


1.16 



Table 4. Estimated mass yields for O, Si and Fe. All values are 
in solar masses. 



Element 


Estimated 
mass 


Type 1A 


Type II 


W7 


WDD2 


12M Q 


2OM 


O 


0.27 


0.14 


0.07 


0.22 


1.48 


Si 


0.032 


0.15 


0.27 


0.099 


0.095 


Fe 


0.033 


0.70 


0.67 


0.059 


0.075 



Next, we check if the total Fe and O masses are consis- 
tent with a Type la SN. According to our volume estimates 
for the Fe-L and Fe-K components (2 and 3) in the previous 
section, their emission comes from volumes of V = 0.85 and 
1.2 xlO 51 m 3 , respectively. Using these volumes along with 
the derived EM and Fe abundance, we calculate an Fe mass of 
0.02 M Q (for / = 0.5) from Fe-L and 0.013 M from Fe-K, 
giving a total Fe mass of 0.033 M Q . Similarly, we use the vol- 
ume estimates made in section 4.2 , the emission measure and 



abundances of components 1 & 2 and 2 & 3 to calculate the 
total O and Si masses respectively. We obtain a total O mass 
of 0.27 M Q and Si mass of 0.032 M©. In table Qwe compare 
the estimated O, Si and Fe masses to model predicted yields for 
type la and type II SNe. The estimated masses are not in exact 
agreement with any of the models, but compare better to the 
type II than type la models. In particular, the Fe mass estimate 
is only a factor 2-3 smaller than predicted by the core collapse 
models, while being more than an order of magnitude smaller 
than can be expected for a type la SN. The mass increases only 
with the square root of the filling factor (or by W), so it is un- 
likely that our simple volume estimates could cause the order 
of magnitude discrepancy between our Fe mass estimate and 
that expected from the Type la models. 
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5. Summary & conclusions 

The high resolution XMM-Newton emission-line spectra of 
SNR N 103B can be fitted well with three plasma components. 
The dominant component (Fe-L) is already in equilibrium and 
can be represented by a kT c = 0.65 keV CIE model. The fact 
that 0.65 keV gas has equilibrated, together with its volume es- 
timate, provides an approximate lower limit to the age of SNR 
N 103B, which is found to be 3000 yrs, twice as old as pre- 
viously estimated. A separate component is needed to explain 
the O-K emission and here we find difficulty to distinguish be- 
tween a relatively cool medium (kT c — 0.2 keV) in equilib- 
rium or hotter (kT c = 0.55 keV) ionising gas. Both scenarios 
produce essentially the same spectrum for O VII. The Fe-K 
emission can be represented by a hot kT c = 3.5 keV NEI com- 
ponent with a relatively low ionisation age of n c i~5.3xl0 16 
m~ 3 s, implying recent shock heating (~ 200 yrs). The de- 
scription of the X-ray emission from SNR N 103B in terms of 
this three-component model is strongly supported by the differ- 
ent morphologies revealed with the Chandra images in three 
corresponding energy bands. The O-K component is concen- 
trated in the northwest, the Fe-L component extends along the 
western rim of the remnant, and the hottest component (Fe-K) 
peaks at a bright knot in the southwest. The O-K and Fe-K 
components lie, at least partially, outside of the Fe-L compo- 
nent, which could imply that the Fe-L component was shocked 
earlier, whereas the O-K and Fe-K regions were shocked only 
recently. However, the different locations and different temper- 
atures of the three components suggest that they are probably 
unrelated dynamically. 

The RGS spectrum unambiguously reveals emission from 
species of O, Ne and Mg. The presence of these lines and 
the abundance measurements they provide suggest that SNR 
N 103B might actually be the result of a type II SN and not a 
type la as previously thought. Also, the total O mass is higher 
than predicted by a type la SNe, while the total Fe mass is much 
lower than would be expected if SNR N 103B was the result of 
a type la SN. 
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